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•rf ; ABSTRACT 

' One third of present-day spirals host optically visible strong bars that drive their dynamical evolution. 

However, the fundamental question of how bars evolve over cosmological times has yet to be resolved, and 
even the frequency of bars at intermediate redshifts remains controversial. We investigate the frequency 
^ \ of bars out to z ~ 1 drawing on a sample of 1590 galaxies from the Galaxy Evolution from Morphology and 

SEDs survey, which provides morphologies from Hubble Space Telescope Advanced Camera for Surveys 
(ACS) two-band images and accurate redshifts from the COMBO-17 survey. We identify spiral galaxies 
using three independent techniques based on the Sersic index, concentration parameter, and rest-frame 
color. We characterize bar and disk features by fitting ellipses to F606W and F850LP images, using the 
two bands to minimize shifts in the rest-frame bandpass. We exclude highly inclined (i > 60°) galaxies 
to ensure reliable morphological classifications and apply different completeness cuts of My < —19.3 and 
—20.6. More than 40% of the bars that we detect have semi major axes a < O'.'S and would be easily 
] missed in earlier surveys without the small point spread function of ACS. The bars that we can reliably 

■ detect are fairly strong (with cUipticities e > 0.4) and have a in the range ~1.2-13 kpc. We find that 

' the optical fraction of such strong bars remains at ~30% ± 6% from the present-day out to look-back 

Tj" , times of 2-6 Gyr (z ^ 0.2-0.7) and 6-8 Gyr {z ~ 0.7-1.0); it certainly shows no sign of a drastic decHne 

' at z > 0.7. Our findings of a large and similar bar fraction at these three epochs favor scenarios in which 

cold gravitationally unstable disks are already in place hy z ^ 1 and where on average bars have a long 
lifetime (well in excess of 2 Gyr). The distributions of structural bar properties in the two slices are, 
I ' however, not statistically identical and therefore allow for the possibility that the bar strengths and sizes 

O ^ may evolve over time. 

' Subject headings: galaxies: evolution — galaxies: general — galaxies: spiral — galaxies: structure 

1. INTRODUCTION gel, & Tcubcn 1997), and starbursts (e.g.. Hunt & Malkan 

T, . . , , • J j-u i 4- 11 u vT, 1999: Jogee, Scoville, & Kenney 2004a), suggests that bars 

It IS widely recognized that stellar bars, either spon- , , . n i • i i • U 

rN , , ,.j 11 • J J J- i -u i J strongly mfiuence their host galaxies. However, the most 

taneously or tidally induced, redistribute mass and an- r i , i • i , . i i ^ ttti ^ 

rrt 1 4- -ij fundamental issues have yet to be resolved. When and 

^ , gular momentum and thereby drive the dynamical and , r ^ a ^ , i 

■ " , 1 f 1 • (- J moo QUI how did bars form.'^ Are bars a recent phenomenon or 

secular evolution 01 galaxies (e.g., Kormendy 1982: Shlos- , , ; , , • i o a i 

„ T? 1 p T3 1 inon T>f ■ p AT innn were they abundant at early cosmic epochs.'^ Are bars 

man, trank, & Begelman 1989; Ptenmger & Norman 1990; , t, ,;, i r 

-c, . ,,. (, o inno J c • j-j- r,nnA\ t long-lived Or do they recurrently dissolvc and re-torm ovcr 

Friedh & Benz 1993; Kormendy & Kcnnicutt 2004 . In rf i , i ; • o tt i . i, i n. , • 

1 , . J r 1 1 -1 1 • u i a Hubble time.'^ How do stellar bars fit withm the hierar- 

the local universe, one-third oi local spiral galaxies host i • i i , • r i n i i ; • i i , 

,. ,1 ■ -1 1 , 1 / 1 -J i 1 onno chical clustering framework of galaxy evolution and relate 

optically visible strong bars (e.g., Eskridge et al. 2002, , , ;• o 

1, fj. T7no 1 c o^ ivj -J • 1 J- to the underlying disk evolution.'^ 

hereafter E02; see also g 3). Mounting evidence, including evolution of a bar over a Hubble time de ends on 

observations of central molecular gas concentrations (e.g., , , . , , , , , /t^t> ;r\ i i ^ 

o 1 + + 1 1nnn^ ^ a ^A f tj \r the host galaxy structure, the dark matter (DM) halo, and 

Sakamoto et al. 1999), velocity fields (e.g., Regan, Vo- .° , at • , • , ■• r .^ ■ ^ 

the environment. Numerical simulations of this complex 
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process make widely different predictions (e.g., Fricdli & 
Benz 1993; Shlosman & Noguchi 1993; El-Zant & Shlos- 
man 2002; Athanassoula 2002; Bournaud & Combes 2002; 
Shen & Sellwood 2004), while the handful of observational 
results on bars at intermediate redshifts are conflicting. 
On the basis of 46 moderately inclined galaxies imaged 
with the Wide Field Planetary Camera2 (WFPC2) in the 
Hubble Deep Field (HDF), Abraham et al. (1999, hereafter 
A99) claim a dramatic decline in the rest-frame optical bar 
fraction (/opt) from -24% at z - 0.2 0.7 to below 5% at 
z > 0.7. On the basis of Near-Infrared Camera and Multi- 
Object Spectrometer (NICMOS) images of 95 galaxies in 
the HDF at z - 0.7-1.1, Sheth et al. (2003, hereafter S03) 
detect four large bars with mean semi major axes a of 
12 kpc (1'.'4), while smaller bars presumably escaped de- 
tection because of the large NICMOS point spread func- 
tion (PSF). S03 point out that their observed bar fraction 
of ^5% for large (12 kpc) bars at z > 0.7 is at least compa- 
rable to the local fraction of similarly large bars. From a 
study based on HubbleSpaceTelescope (HST) Advanced 
Camera for Surveys (ACS) F814W images of 186 galax- 
ies over a 3'.9 x 4'.2 area in the Tadpole field, Elmegreen, 
Elmegreen, & Hirst (2004, hereafter E04) report a con- 
stant optical bar fraction of ~20%-40% at z ~ 0-1.1. This 
study is limited by the large (0.1-0.4) errors of the photo- 
metric redshifts (Benitez et al. 2004), and the small sam- 
ple size precludes absolute magnitude completeness cuts. 
Furthermore, with only a single ACS filter, the rest-frame 
bandpass of the observations shifts by more than a factor 
of 2 over z ~ 0-1.1. 

In this Letter (see also Jogce et al. 2004b), we present 
the first results of an extensive study of bars at z ~ 0.2- 
1.0, based on two-band ACS images covering 14' x 14' 
(~25%) of the Galaxy Evolution from Morphology and 
SEDs (GEMS) survey. The area and the sample size of 
1590 galaxies provide at least an order of magnitude better 
number statistics than earlier studies. We quantitatively 
identify bars using ellipse fits and minimize the effects of 
bandpass shifts by using both F850LP and F606W images 
(§ 2.2). Using highly accurate redshifts (§ 2.1), we com- 
pare the bar fractions in two redshift slices after applying 
completeness criteria. We show that the optical fraction of 
strong (ellipticities e > 0.4) bars is remarkably constant (~ 
30%) from the present-day out to look-back times (Tback) 
of 2-6 Gyr (z ~ 0.2-0.7)^^ ^nd 6-8 Gyr (z ~ 0.7-1.0). 

2. OBSERVATIONS, SAMPLE, AND METHODOLOGY 

2.1. Observations and Sample Selection 

GEMS is a two-band (F606W and F850LP) HST ACS 
imaging survey (Rix et al. 2004) of an 800 arcmin^ (~ 
28' X 28') field centered on the Chandra Deep Field-South. 
GEMS consists of 78 one-orbit- long ACS pointings in each 
filter and reaches a limiting 5 a depth for point sources of 
28.3 and 27.1 AB mag in F606W and F850LP, respec- 
tively. GEMS provides high-resolution (— 0'.'5 or 360 pc 
at z — 0.7) ACS images for —8300 galaxies in the red- 
shift range z — 0.2 1.1, where accurate redshifts [5z/(l + 
z) ~ 0.02 down to i?vega = 24] and spectral energy distri- 
butions (SEDs) based on 17 filters exist from the COMBO- 
17 project (Wolf et al. 2004). For this Letter, we analyze 
only ~25% of the GEMS field as this area (14' x 14') is 

We assume in this paper a flat cosmology with CIm = 1 — f^A = 0.3 



already 30 times that of the HDF and yields good num- 
ber statistics and robust results on the bar fraction (§ 3). 
It provides a sample that consists of 1590 galaxies in the 
range z ~ 0.2-1.0 and i?vcga < 24. In a future paper (S. 
Jogee et al. in preparation, hereafter Paper II), we will use 
the entire GEMS sample to compare how bar properties 
evolve over 8 Gyr at 1 Gyr intervals. 



2.2. Characterizing bars and /opt out to z 1.0 

Table 1 illustrates the two methods that we use for 
assessing bar properties in two redshift slices at 0.25 < 
z < 0.70 and 0.7 < z < 1.0. The first method (referred to 
as Method I in Table 1) is to identify bars in the reddest 
filter (F850LP) at all z in order to minimize extinction 
and better trace old stars. However, with this method, 
the rest-frame bandpass shifts significantly, from I to B 
across the redshift range 0.2-1.0. The second complemen- 
tary method is to trace bars in both F606W and F850LP 
images such that the rest-frame band remains relatively 
constant, between B and V, out to z ~ 1.0. 

We identify bars and other galactic components in 
F606W and F850LP images via the widely used (e.g., 
Wozniak et al. 1995; Jogee et al. 1999, 2002; Knapen et 
al. 2000) procedure of fitting ellipses using the standard 
IRAF "ellipse" routine. We developed a wrapper that au- 
tomatically runs "ellipse" for a range of different initial pa- 
rameters, performing up to 200 fits per galaxy. We success- 
fully fitted ellipses to 90% of the 1590 galaxies, while the 
10% failure cases included mostly very disturbed systems 
where no centering could be performed and some extended 
low surface brightness systems. For all fitted galaxies, we 
inspected the image (Fig. la), the fitted ellipses overlaid 
on the images (Fig. lb), and the radial profiles (Fig. Ic) 
of intensity, ellipticity (e), and position angle (P-A.) in or- 
der to confirm that the best fit is reliable. This inspection 
was aided by a visualization tool that we developed. We 
identify a bar as such if the fitted ellipses and radial pro- 
files show the following characteristic bar signature (e.g, 
Wozniak et al. 1995) illustrated in Fig. 1. (i) The ellip- 
ticity (e) must rise to a global maximum Cmax, which we 
require to be above 0.25, as well as above that of the outer 
disk, while the P. A. has a plateau (within ± 20°) along 
the bar. (ii) Beyond the bar end, as the bar-to-disk tran- 
sition occurs, e must drop by > 0.1, while the P. A. usually 
changes by > 10°. From the profiles, we identify e, P.A., 
and semi major axes a of both the bar and the outer disk. 
We quantify the bar strength using emax, which correlates 
locally with other measures of bar strength, such as the 
gravitational bar torque (Laurikainen, Salo, & Rautiainen 
2002). 

The sample of 1430 galaxies with successful ellipse fits 
includes galaxies of different morphological types (disks 
and spheroids), inclinations, and absolute magnitudes. We 
apply two cuts at My < —19.3 and —20.6. The first cut 
gives us a sizeable sample of galaxies with a range of ab- 
solute magnitudes (—19.3 to —23.8) similar to that of the 
Ohio State University (OSU) survey, which is used to de- 
fine the local bar fraction (E02). However, ii'-corrections 
based on local Scd templates (Coleman, Wu, & Weedman 
1980) suggest that we are only complete out to z ~ 0.8 for 

and Ho =70 km s"^ Mpc~^. 
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the first cut. The second more stringent cut at —20.6 en- 
sures completeness out to z ^ 1.0, but it reduces the sam- 
ple drastically (see Table 1). In addition, to ensure reliable 
bar detection, we use the disk inclination i from ellipse fits 
(§ 2.2) to exclude highly inclined {i > 60°) galaxies. 

The optical bar fraction /opt is defined as 
(-^Vbar/-/Vsp/disk) where A?sp/disk is the number of spiral or 
disk galaxies, and A^bar is the number of such systems host- 
ing bars. We identify spiral/disk galaxies using three in- 
dependent techniques (Table 1). We first use the criterion 
n < 2.5, where n is the index of single-component Sersic 
models fitted to GEMS galaxies (B. Haussler et al. 2004, 
in preparation) with the GALFIT (Peng et al. 2002) pack- 
age. Our choice of n < 2.5 is dictated by the fact that a 
study of GEMS galaxies at z - 0.7 (Bell et al. 2004), as 
well as tests in which we insert artificial galaxies in the 
GEMS fields, suggest that a Sersic cut of n < 2.5 picks 
up disk-dominated systems. The second technique uses a 
cut C < 3.4, where C is the CAS (Consehce et al. 2000) 
concentration index. Third, we use rest-frame U — V color 
cuts in the range 0.8-1.2 to broadly separate spiral galax- 
ies from red E/SOs, based on local SED templates and the 
observed red sequence at ^ ^ 0.7 (Bell et al. 2004). 

3. RESULTS AND DISCUSSION 

The bars that we identify primarily have ellipticities 

e > 0.4 and semi major axes a in the range 0'.'15-2'.'2 and 
1.2-13 kpc (Fig. 2). Our experiment of artificially red- 
shifting B-band images of a subset of OSU galaxies out to 
z ^ 1 shows that it is difficult to unambiguously identify 
weaker (e < 0.3) bars, and we limit the discussion in this 
Letter to strong (e > 0.4) bars. Table 1 shows the opti- 
cal bar fraction /opt of such bars in the two redshift slices 
(0.25 < z <0.70 and 0.7< z <1.0) derived in the reddest 
filter and in a relatively fixed rest-frame band (§ 2.2). Re- 
sults for My < —19.3 are shown, based on 627 galaxies 
out of which we identify 258 moderately inclined {i < 60°) 
spirals that host ~80 bars. The consistency in the six en- 
tries attests to the robustness of the results and shows that 
the fraction of optically visible bars remains in the range 
23%-36% or at ~30% ± 6% in both slices. Incomplete- 
ness effects (§ 2.2) do not bias the results since the cut at 
—20.6 gives similar bar fractions, shown in brackets (Table 
1). The bar fraction is slightly higher in the rest-frame / 
band than B band, possibly indicative of dust and star 
formation masking bars at bluer wavelengths. Our find- 
ings of a fairly constant /opt are consistent with E04 and 
do not show the dramatic decline in /opt reported by A99. 

More than 40% of the bars that we detect have semi 
major axes a < 0'.'5 (Fig. 2), and many of these small- 
est bars may have been missed in earlier WFPC2 (e.g., 
A99) and NICMOS studies that did not benefit from the 
small (0'.'05) ACS pixels, and the resulting narrow PSF. In 
addition to the wider WFPC2 PSF, cosmic variance, low 
number statistics, and methodology may have led to the 
lower /opt reported by A99, but we cannot address this 
issue further here as the coordinates of galaxies in that 
study have not been published. 

How do our results compare to /opt for correspondingly 
strong bars in the local universe? There are as yet no 
statistics published on ^; bars based on a large volume- 
limited sample such as the Sloan Digitized Sky Survey. 



As the next best alternative, we turn to the OSU sam- 
ple (E02). For i < 60° spirals, we find /opt 37% for 
strong bars classified according to the visual RC3 'SB' bar 
class, and /opt ~ 33% for strong bars classified accord- 
ing to e > 0.4, where e is based on ellipse fits to OSU 
i?-band images. Thus, it appears that the optical fraction 
of strong (e > 0.4) bars remains remarkably similar at 
^30%-37% from the present-day out to look-back times of 
2-6 Gyr (z ~ 0.2-0.7) and 6-8 Gyr (z ~ 0.7-1.0). Our 
findings have several implications for disk and bar evo- 
lution. (1) The abundance of strong bars at early times 
implies that dynamically cold disks that can form large- 
scale stellar bars are already in place by z ^ 1. They also 
suggest that highly triaxial, centrally concentrated DM 
halos, which tend to destabilize the bar (El-Zant & Shlos- 
man 2002), may not be prevalent in galaxies at z ~ 0-1. 
(3) The remarkably similar fraction of strong bars at look- 
back times of 6-8 Gyr, 2 6 Gyr, and in the present day 
supports the view that large-scale stellar bars are long- 
lived (e.g., Athanassoula 2002; Shen & Sellwood 2004; 
Martinez- Valpuesta et al. 2004), with a lifetime well above 
2 Gyr. The alternative option, statistically allowed by the 
data, is that the combined destruction and reformation 
timescale of bars is on average well below 2 Gyr. How- 
ever, this is highly implausible because the destruction of 
a bar leaves behind a dynamically hot disk that cannot 
reform a bar unless it is substantially cooled via processes 
that take at least several Gyr, such as the accretion of 
large amounts of cold gas. (3) We note that a similar frac- 
tion of bars at different epochs does not exclude the pos- 
sibility that the bar strength (ellipticity) and its size can 
evolve in time because of intrinsic factors and concurrent 
changes in the surrounding disk, bulge and halo potentials 
(e.g., Athanassoula 2002; Martinez- Valpuesta & Shlosman 
2004). Kolmogorov-Smirnov tests on the distributions of e 
and a for different magnitude cuts yield primarily P in the 
range 0.2-0.5, where P is the significance level for the null 
hypothesis that the two data sets are drawn from the same 
distribution. Such values of P hint at evolutionary effects, 
but a larger sample is needed to draw definite conclusions. 
While this Letter focuses on the optical bar fraction, in 
Paper II, we will invoke the full GEMS sample of 8300 
galaxies to compare bar and host galaxy properties (e.g., 
disk scale lengths, masses, and B/D ratios) over the last 
8 Gyr at 1 Gyr intervals. 
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Table 1 

Optical fraction /opt of strong (e >0.4) bars at 
0.25 <z< 0.70 (Tback ~ 2-6 Gyr) and 0.7 < z < 1.0 (Tback ~ 6-8 Gyr) 



Redshift 


ATgal 


Technique to 


-^sp/disk 


Filter 


Rcsl-rrame 


Optical bar 


range 




identify 


to trace 


band 


fraction 




disks/spirals 




bars 




/opt 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



Method I: Using the reddest ACS filter (F850LP) to identify bars 

0.25-0.70 384 (146) Sersic n 148 (39) F850LP / to y 36% (33%) 

0.70-1.0 243 (13,5) Sersic n 110 (49) F850LP V to B 24% (27%) 

0.25-0.70 384 (146) CAS C 170 (51) F850LP / to y 27% (27%) 

0.70-1.0 243 (135) CAS C 147 (73) F850LP y to B 23% (23%) 

0.25-0.70 384 (146) U -V 175 (51) F850LP / to y 29%(28%) 

0.70-1.0 243 (135) U -V 139 (69) F850LP V to B 29% (26%) 



Method II: 


Using F606W and F850LP to 


1 identify 


bars at an 


approximately fixed rest-frame B/V band 


0.25-0.70 


384 (146) 


Sersic n 


148 (39) 


F606W 


y to B 


26% (30%) 


0.70-1.0 


243 (135) 


Sersic n 


110 (49) 


F850LP 


y to B 


24% (27%) 


0.25-0.70 


384 (146) 


CAS C 


170 (51) 


F606W 


y to B 


24% (28%) 


0.70-1.0 


243 (135) 


CAS C 


147 (73) 


F850LP 


y to B 


23% (23%) 


0.25-0.70 


384 (146) 


U - V 


175 (51) 


F606W 


y to B 


29% (23%) 


0.70-1.0 


243 (135) 


U - V 


139 (69) 


F850LP 


y to B 


29% (26%) 



Note. — Col. (1): The redshift range. Col. (2): The no of galaxies in this range with My < —19.3 and 

— 20.6 (shown in brackets) to which ellipses were fitted. Col. (3): The technique used to identify disk/spiral 
galaxies from the sample in (2). We use cuts of Sersic index n < 2.5, CAS concentration index C < 3.4, and 
rest-frame U — V < 0.8—1.2. Col. (4): A'^gp/diski of moderately inclined {i < 60°) spiral/disk galaxies for 

the two magnitude cuts in (2). Col. (5), (6): The filter and rest-frame band in which bars are traced; Col. (7): 
/opt, the optical fraction of strong bars, is the fraction (A''bar/-^sp/disk) of moderately inclined spirals hosting 
bars with e > 0.4. Values for the two magnitude cuts of —19.3 and —20.6 are shown. 
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Fig. 1. — Characterization of bars out to ^ ~ 1.0: The GEMS image of a 2 ~ 0.5 galaxy with a bar, prominent spiral arms, and a disk is 
shown without (a) and with (b) an overlay of the fitted ellipses, (c) In the resulting radial plots of the surface brightness, ellipticity e, and 
P. A., the bar causes e to rise smoothly to a global maximum, while the P.A. has a plateau. Beyond the bar end (a ~ 0'.'36), the spiral arms 
lead to a twist in P.A. and varviug c ijoforc the the disk dominates. 
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Fig. 2. — Compajrison of bars out to look-back times of 8 Gyr: The bar ellipticity e (left) and semi major axis a in " (middle) and 
kpc (right) are shown for bright (My ^ —19.3), moderately inclined (i < 60°) galctxies at ^ ~ 0.2-0.7 {T^^sA ~ 2-6 Gyr) and z ~ 0.7-1.0 
{^back~ 6~8 Gyr). The bars identified are primarily strong, with e > 0.4. A large fraction have a < 0'.'5 and their detection is aided by the 
narrow ACS PSF. 



